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The paper contains further development of the idea of field quantization introduced in M. Cza- 
chor, J. Phys. A: Math. Gen. 33 (2000) 8081-8103. The formalism is extended to the relativistic 
domain. The link to the standard theory is obtained via a thermodynamic limit. Unitary repre- 
sentations of the Poincare group at the level of fields and states are explicitly given. Non-canonical 
multi-photon and coherent states are introduced. In the thermodynamic limit the statistics of pho- 
tons in a coherent state is Poissonian. The S matrix of radiation fields produced by a classical 
current is given by a non-canonical coherent-state displacement operator, a fact automatically elim- 
inating the infrared catastrophe. Field operators are shown to be operators and not operator- valued 
distributions, and can be multiplied at the same point in configuration space. An exactly solvable 
example is used to compare predictions of the standard theory with those of non-canonical quan- 
tum optics, and explicitly shows the mechanism of automatic ultraviolet regularization occuring 
in the non-canonical theory. Similar conclusions are obtained in perturbation theory, where one 
finds the standard Feynman diagrams, but the Feynman rules are modified. A comparison with the 
Dicke-Hepp-Lieb model allows to identify the physical structure behind the non-canonical algebra 
as corresponding to an ensemble of indefinite- frequency oscillators with constant density N/V. 

PACS: ll.lO.-z, 04.60. Ds, 98.80.Es 



I. INTRODUCTION 

The idea of "non-canonical quantization" of electromagnetic fields was introduced in Q. One begins with the 
observation that even in nonrelativistic quantum mechanics it is natural to treat the frequency uj characterizing a 
harmonic oscillator as an eigenvalue and not a parameter. A replacement of uj by an operator lj leads to an indefinite- 
frequency oscillator with altered (non-)canonical commutation relations, and forms a natural departure point for a 
new version of the 'old fashioned' field quantization. An analysis of physical structures associated with such indefinite- 
frequency operators was discussed in detail in |^ . In the present paper we show that the non-canonically quantized 
electromagnetic field may be regarded as an ensemble of indefinite-frequency oscillators of constant density. 

We extend the formalism to the relativistic domain. Representations of both the non-canonical commutation re- 
lations and the Poincare group are explicitly constructed. Multi-photon and coherent states are defined. Their 
properties make them similar to those from the canonical formalism if one performs an appropriately defined thermo- 
dynamic limit N ^ oo. Vacuum states transform as a massless scalar field. Radiation fields produced by a classical 
current lead to the correct form of the S matrix. The S matrix is proportional to a non-canonical coherent-state 
displacement operator, a fact eliminating the infrared catastrophe. 

Commutation relations for vector potentials are found showing certain deviations from locality due to nontrivial 
structure of vacua but simultaneously allow for multiplication of field operators at the same point in configuration 
space, a first hint suggesting the ultraviolet finiteness of the theory. Two other arguments in favor of the ultraviolet 
finiteness are based on analysis of perturbation theory and the spontaneous emission survival amplitude. 

An example of the spontaneous-emission survival amplitude is exactly solvable in both canonical and non-canonical 
theories and can be employed to illustrate the links between the two approaches. Of particular interest is the analogy 
between our method of quantization and the Hepp-Lieb treatment of the Dicke model. We show that that the Dicke- 
Hepp-Lieb Hamiltonian has the structure of the same type as our multi-oscillator non-canonical Hamiltonian. This 
leads to the conclusion that the formal requirement that the RHS of non-CCR commutators satisfies the resolution 
of unity is equivalent to the physical requirement that the electromagnetic field is an ensemble of indefinite-frequency 
harmonic oscillators of constant density N/V. 

The fact that the RHS of commutation relations is modified by the presence of an element from the center of 
the algebra makes our "non-canonical" fields analogous to the so-called generalized free fields [^,||. The essential 
difference between the two approaches is that, first of all, our fields do interact with charges and we do not demand 
the Poincare invariance of field commutators but employ the algebra which is only covariant. As a result we loose 
Poincare invariance of vacuum and locality. This leads to analogies with nonlocal quantum field theories discussed 
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in . The loss of locality implies also interesting analogies between our approach and quantum field theory in non- 
commutative space |^-|^. Still, there are also differences between what we propose and these approaches. For example, 
the formalism of is based on nonlocal distributions and regularizations of Hamiltonians. Similarly to our approach 
vacua in such theories are Poincare non- invariant, a fact used to circumvent limitations imposed by the Haag theorem 
1^. However, in our approach the nonlocal distributions occur at the level of amplitudes or averages, and not at the 
algebra of field operators. Analogously, perturbative expansions of amplitudes are automatically regularized in spite 
of the fact that there is no regularization at the level of operators. 



II. NOTATION 



In order to control covariance properties of fields in generalized frameworks it is best to work in a manifestly 
covariant formalism. The most convenient is the one based on spinors and passive unitary transformations. 



A. Spinor convention and fields 



We take c = 1 and h = I. The index notation we use in the paper is consistent with the Penrose- Rindler spinor 
and world- tensor convention ^J. The electromagnetic field-tensor and its dual are 

El E2 

p I -El -02 I , s 

- I -Bi I 

— i?3 —B2 Bi 

~Bi -B2 -Bi 
p _ I -Bi -£'3 E2 I 
^"fc - \ B2 E:i -El I 

i?3 — E2 El 

Self-dual and anti-self-dual parts of Fab are related to the electromagnetic spinor by 

-Fi F2 F3 
l^r. \ . - { -Fi iFs -iF2 



'Fab = -\^Fab-i*Fabj =eABVA'B' - \ _ _-p^ q -p^ \ (3) 



-F3 iF2 -iFi 



-F 



^ (^Fab + i*Fab^ = eA'B'VAB (4) 

where F — {E + iB)/2 is the Riemann-Silberstein vector Denote k ■ x — kax"" . The electromagnetic spinor has 
the following Fourier representation [pT 12 



^ab{x)^ j dT{k)^A{k)^B{k)[f{k,-)e-'^-^ + f{k,+)e'^--') (5) 

^A'B' {x) = J dr{k)nA' (fc)^B' (k) (/(fc, -H)e-*^-^ -t- Jii^e^''-^) (6) 

where the spinor field 7r^(fc) is related to the future-pointing 4- momentum by 

k" = 7r^(fc)7f^'(fe) = (fco,fc) = (|fc|,fc) (7) 

and the invariant measure on the light-cone is dr(k) — [(27r)'^2fco] ^cPk. Anti-self-dual and self-dual parts of the field 
tensor are 

-Fabix) = J dr(fe)e^,B,7r^(fc)7rB(fc)(/(fc,-)e-'''^--^ +7(M:)e*^) (8) 

= I drik)eabik) (/(fc, -)e-*- + 7(fc;T)e*-) (9) 

+Fab{x) = I dT{k)eAB7TA'{k)7tB'{k)(jik;^e'''-^ + /(fe,+)e-*-) (10) 

drik)eabik) (7(fc~)e*- + /(fc, +)e-*-) (11) 
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The latter form is used by the Bialynicki-Birulas in jl^ . 

The sign in the ampUtude /(fc, ±) corresponds to the value of helicity of positive-frequency fields. 
The four- vector potential Aa{x) is related to the electromagnetic spinor by 

Vxy{x) ^V(^x^'Ay)y'{x) (12) 
In the Lorenz gauge V^Aa = we do not have to symmetrize the unprimed indices and 

ipXY{x)^Vx^'AYY'ix). (13) 

One of the possible Lorenz gauges is 

Aaix) = z y dr(fc)|^ma(fc)(/(fc,+)e-*^ -7(fc~)e*-) - Tfi.lfe) (7(fc;T)e*" - /(fe,-)e-*-)^ (14) 

where lUj^tt"^ = 1, i.c lua = LUA{k) is a spin- frame partner of Tr^{k). In ( p^ wc have introduced the null vectors 

ma{k) ^ uJA{k)nA'{k) (15) 

ffiaik) TrA{k)LJA'{k) (16) 

which, together with 

ka = nA{k)7fA'{k) (17) 

UJa{k) = LJA{k)u!A'{k) (18) 

form a null tetrad |^ . 

A change of gauge is in the Fourier domain represented by a shift by a multiple of k"^. The form (^4|) shows that 
gauge freedom is related to the nonuniqueness of u}A{k) which can be shited by a multiple of TTA{k). 

B. Momentum representation 

Consider the momentum-space basis normalized by 

{p\p') - (27r)32po^''Hp-p') - Srip,p'), Po > 0. (19) 

The identity operator in momentum space is / dr{p)\p){p\. We can use the following explicit realization of |p) = \fp) 
in terms of distributions 

fp{k) = (2^)32po<5('^(p -k)= 6r{p, k) (20) 

Since 

J dT{k)Fik)Sr{p,k) = J d^k5^^\p - k)F{k) = F{p) (21) 
the Fourier transform of fp{k) is 

Uix) = J dr(fc)/p(fc)e'*^ = e-'P-^ (22) 

If 1 is the identity operator occuring at the right-hand-side of CCR [a^, aj,] = (5ss'l, we denote 

/fc = |fc)(fe|0l, (23) 
/ = / dr{k)\k){k\^l. (24) 
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C. Multi-particle conventions 



Let A be an operator A : H ^ H where H is a one-particle Hilbert space. The muhi-particle Hilbert space 

2i = ®;f=iCW (25) 

is the Hilbert space of states corresponding to an indefinite number of bosonic particles; ®"7i stands for a space of 
symmetric states in H . ^ (^H . We introduce the following notation for operators defined at the multi-particle level: 

n 

©a„ A = aiA ® a2{A I + I A) ® a3{A IS) I I + I A I + I I A) ® . . . (26) 

Here ®a„ A : H ^ an are real or complex parameters, and / is the identity operator in H. 
The following properties follow directly from the definition 

[©„„ A, ®/3„B] = ®a„/3,. [A, B] (27) 

ge„„A ^ e^^^ e""^ e""-^^ (28) 

n 

e®^-^ ©ft. Be-®'"" = (Dp^e^Be-^ (29) 
Identity operators in H and H are related by 

/ = ®i/. (30) 

We will often use the operator 

Lk = ®±Ik. (31) 



III. POINCARE TRANSFORMATIONS OF CLASSICAL ELECTROMAGNETIC FIELDS 

Denote, respectively, by A and y the S'L(2, C) and 4-translation parts of a Poincare transformation iQ (A, y). 
The spinor representation of the Poincare group acts in the space of anti-self-dual electromagnetic fields in 4-position 
representation as follows: 

-Fat{x)^ {TA^y-F)Jx) (32) 
= A/Ab'^-/',rf(A-i(x-y)) (33) 

= J dr(fe)£A'B'AA^^c(fc)AB^7rz3(fc)(/(fc,-)e-*^~'(--^)+7(fc:T)e*''"'("-^)) (34) 
= J dr(fe)£A'B'AA^^c(A-'fc)AB^^z3(A-ifc)(/(A-ife,-)e-*("-^)+/(A-ifc,+)e*(^-^)) (35) 

where A~^k is the spacelike part of A^^a^'fcb- The transformed field 

(A^)^(fc) = A^^^c(A-^fc) (36) 

satisfies 

fc" ^ Tr^{k)T:^'{k) = (A7r)'4(fc)]U^'(fc) (37) 

Now, if uJA{k) is a spin-frame partner of 7r^(A;), i.e. aj^(A;)7r'^(A;) — 1, one can write 

TT^{k) ^ {An)^{k)QA'{k)A^'^'{k) (38) 

which shows that n^{k) and (A7r)^(fc) = A^'-^7r(7(A~-'-fe) are proportional to each other, the proportionality factor 
being 
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X{A,k) =u;A'ik)ATr^'ik). (39) 
The form showed that the gauge freedom is related to shifts 

uJA{k) y—f tJA{k) + scalar x TTA{k) (40) 
which do not affect A(A, k) making it independent of gauge. Using again 

r = |A(A,fc)|27r'^(fe)7f-^'(fc) = |A(A,fe)|2fc'» (41) 
one concludes that A(A, k) is a phase factor 

A(A, fe) = e*®('^^'^) (42) 

and we find 

{TA,y~F)Jx) = I dT{k)eA'B'^A{k)^B{k)e-^''"^^^''^ (/(A'^fe, -)e-*(--^) + /(A'^fe, +)e*(--^)) (43) 

dr(fe)eafc(fc)((rA,^/)(fc,-)e-*- + (TA,,/)(fc,+)e*-) (45) 

We have obtained therefore the passive transformation of the classical wave function 

/(fc, ±) ^ (T^,yf){K ±) = e±2»^(^^'=)e*V(A-'fc, ±)- (46) 

which is simply the unitary zero-mass spin-1 representation the Poincare group. The above derivation clearly shows 
that the rule ( ^6| ) is obtained without any particular assumption about the choice of /(fc,±). In particular, the 
derivation remains valid even if one replaces functions / by operators, independently of their algebraic properties. The 
above 'passive' viewpoint on the structure of unitary representations is particularly useful if one aims at generalizations 
of CCR. The passive derivation of all the non-tachyonic unitary representations of the Poincare group can be found 
in Hi]. 

IV. NON-CANONICAL QUANTIZATION 

We follow the strategy described in |Q] and Let be canonical annihilation operators satisfying CCR [as, a\,] — 
(5ss'l. Define the 1-oscillator non- canonical creation and annihilation operators [p5 | 

«(/)^ =E / rfr(fc)/(fc,s)|fc)(fc|®at (47) 

s •' 

= J2 f dr{k)f{k,s)a{k,sy (48) 
aif) = E / dr{k)J{k;sj\k) {k\<g>as (49) 

= E y dr(fc)7(M)«(fc, s) (50) 



satisfying the non-CCR algebra 



[a(fc,s),a(fc',s')^] ^Sss'Sr{k,k')\k){k\(g)l (51) 
^6ss'5r{k,k')h. (52) 



Taking, in particular, fp^r{k, s) — Srs^rip, k) one finds 

a{fp,r) = \p){p\<»ar=a{p,r). (53) 
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The one-oscillator quantization is 

^ Fabix) = EA'B'^ABix) (54) 

•.r(.)e..,..(.)..(fc)(a(/c,-)e-- + a(.,+)te-) (55) 

dT{k)eab{k) (a(fe, -)e-*^-" + a(fe, +)^e*") (56) 

Spinor transformations of ~ Fab{x) lead to the passive transformation 

a{k, ±) ^ (TA,aa)(fc, ±) = e±2^®(^'*=)e''=-''a(A-ifc, ±). (57) 

The quantization procedure is gauge independent since we work at the gauge- independent level of ~ Fab{x). 
Multi-oscillator fields are defined in terms of 

a(p,s) ®_i_a{p,s) (58) 

and 

= E / '^r(fc)/(fc, s)a{k, s)t (59) 

a(/) = E / c?r(fc)7(fc:7)a(fe, s). (60) 

The fact that the coefficients are found in the multi-oscillator definition may appear awkward. At least three 
different formal arguments for such a choice of the multi-particle extension of field operators were given in j^]. Below, 
in Sec. X, we will show that this very special choice of the non-CCR representation corresponds physically to an 
ensemble of oscillators uniformly distributed in space. 
The non-CCR algebra is 

[a(/), «(5)^] = E / dT{k)J{k;s)gik, s)Lk (61) 

The right-hand-side of the above formula is in the center of the non-CCR algebra, i.e. 

[[aif),aig)%aih)]^0 (62) 

[[a{f).aig)%a{hy]=0 (63) 

Useful is also the formula 

[a{fp,r),a{fp,y)^] - [a(p,r),a(p',r')1 (64) 

^Srr'Sr{p,p')Lp. (65) 

The presence of /p at the right-hand-sides of non-CCR will infiucnce orthogonality properties of multi-photon states, 
as we shall see later. 

At the multi-oscillator level the electromagnetic field tensor operator is 



F^,{x) = J dV{k)eab{k) (a{k, -)e-*^ + a(fe, +)te*^) (66) 

+ j dT{k)eab{k) (aik, -)te*- + a{k, +)e-*") . (67) 
The four-potential operator is (in our choice of gauge) 

^,.,../.r,.,(,M.,(.(.,.)e---.(.,-,.-'i...(.)(.,.-,e---.^ m 

It is well known that field "operators" of the standard theory (let us denote them by Aa{x)) are in fact operator- valued 
distributions. As a consequence the operator products of the form Aa{x)Ab{y) are ill defined and lead to ultraviolet- 
divergent expressions ii x = y. The techniques of dealing with ultraviolet divergences are based on appropriate 
regularizations of products of distributions taken at "diagonals" in configuration space |l^,|l^. In Sec. VI w shall see 
that the non-canonically quantized A^{x) is an operator and there is no difficulty with A^Xx)A^^{x). 
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V. ACTION OF THE POINCARE GROUP ON FIELD OPERATORS 



We are interested in finding the representation of the group in terms of unitary similarity transformations, i.e. 

a(fc, ±) ^ e±2«(A^'=)e*fc-i^a(A-ifc, ±) = {/i_^a(fe, i)!/^,, (69) 

It is sufficient to find an appropriate representation at the one-oscillator level. Indeed, assume we have found Uh,y 
satisfying 

g±2ze(A,fc)g./c.y^(^-i^^ ±) = Ul ya{k, ±)C/a,,. (70) 

Then 

OO 

H^y = yA,y®- - -®UA,y^ ■ (71) 



A. Four-translations 

The definition of four momentum for a single harmonic oscillator is 



where 



One immediately verifies that 



P^^ j dV{k)ka\k){k\®h (72) 
h=]^'^{alas + asa^)j ^^hs. (73) 



a(fc,s)e-'-f'-^ = a(fc,s)e-"-'= (74) 
a(fc, s)^e-''P-^ = a(fc, s)^e"-'= (75) 



implying 



Ui^y - e^y P. (76) 

Consequently, the generator of four-translations corresponding to C/j^ y = e*^ — is = (BiPa and 

e -■"a(fc, s)^e-'^-" = a{k, s)te"-'= (77) 

e'-e-^a(fc, s)e-*£-^ = a(fc, s)e-"-^ (78) 
The x-dependence of field operators can be introduced via P: 

EM = e^^-^F^^- (79) 



B. Rotations and boosts 

To find an analogous representation of 

a{k, ±) ^ e=^2''®(^''=)a(A-^fc, ±) = Hl^oaik, ±)Ua,o (80) 

we define 

t/A,o = exp ( ^ 2is / dT{k)e{A,k)\k){k\(E)hs^(^J2 I dT{p)\p,r){K-'^p,r\ ® I 

s r •' 

Finally the transformations of the field tensor are 

uloPabWAfi = A/Ab'^F,,(A-ix) (82) 

Ul.,yF,,ix)U^^y ^ F,,ix - y) (83) 

The zero-energy part of P_ can be removed by a unitary transformation leading to a vacuum picture dynamics (cf. 
0). We will describe this in more detail after having discussed the properties of non-canonical states. 
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VI. STATES AND THEIR POINCARE TRANSFORMATIONS 



It is clear that in order to control transformation properties of states it is sufficient to discuss single-oscillator 
representations. We shall start with single-oscillator states and then extend them to many oscillators. 

A. Representation in the one-oscillator sector 

The one-oscillator Hilbert space consists of functions / satisfying 

^ /dr(fc)|/(fc,n+,n_)|2 <oo. (84) 

n-j- ,n- —0 

We will write them in the Dirac notation as 

I/) = E / ^r(fc)/(fc, n+, n_)|fc, n+, (85) 

The representation of the Poincare group is 
|/)^C/a,,|/> -C/i,,C/a,o|/) 

= Z1 / '^r(fc)/(A-^fc,n+,7i_)e2*("+-"-)®(^''^)e*^("++"-+i/2)|fc,n+,7i_). (86) 

The latter formula can be written as 

/(fc, n+, n_) ^ UA,vfik, n+,n_) = e*'=-^("++"-+i/2)e2''("+-"-)e(A.fc) /(A'lfe, n+, n_) (87) 

or 

UaJI) = \UA,yf). (88) 
The form (p7|) is very similar to the zero- mass spin-1 representation (Eq), the difference being in the multiplier 



n+ + n_ + 1/2. One can check by a straightforward calculation that ( |87| ) defines a representation of the group. 

B. Generators and vacuum picture 

Denote by Ka and Lab + Sab the generators of 4-translations and SL{2, C) of the standard zero-mass spin-1 unitary 
representation of the Poincare group. Lab denotes the orbital part of the generator. The generators of ( p7| ) are then 

^Ka®h (89) 
^L''^®l + S''\s®hs. (90) 

S""^ ss' are matrix elements of S"'^ (which is a diagonal p-dependent matrix). Denote by Sab the generators of the 
(1/2, 1/2) spinor representation of SL{2,C), i.e. A = exp (i£^°-''Sab/2). The generators of the unitary representation 
are defined by 

_9 

dy 
d 

Jabf{k,n+,n-) = ~i—^UK,yf{k,n+,n-)\^^^^^ (92) 

In what follows we will work in a "vacuum picture" , i.e with unitary transformations 

/(fc,7i+,n_) ^ VA,yfik,72+,n_) = e^("++"-)'=-2'e2*("+'"-)®(^''=)/(A-^fc,n+,n_). (93) 

The transition 

Ua^V ^ VA^y = W^UA,y (94) 

is performed by means of the unitary transformation which commutes with non-CCR creation and annihilation 
operators. 

Let us stress that the fact that we "remove" the zero-energy parts from generators does not mean that energy of 
vacuum is zero. The vacuum picture is in a sense a choice of representation co-moving with vacuum. 



Paf{k,n+,n-) = -i — UA.yf{k,n+,n-)\^^y^^ (91) 
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C. Vacuum states 

Vacuum states are all the states which are annihilated by all annihilation operators. At the one-oscillator level 
these are states of the form 

\0) = j dr{k)0{k)\k, 0,0). (95) 

Even in the vacuum picture the vacuum states are not Poincare invariant since 

VA,yO{k) = O(A-^fe) (96) 

which means they transform as a 4-translation-invariant scalar field. We will often meet the expression Z{k) = |0(A;)|^ 
describing the probability density of the "zero modes" . 

D. Coherent states 

An analogue of the standard coherent (or "semiclassical" ) state is at the 1-oscillator level 

|0„)= J dr{k)0{k)\k)\a{k,+),a{k,-)} (97) 

where 

as|a(fe, +), a{k, -)) = a(fe, s)\a{k, +), a{k, -)) (98) 
Its explicit form in the basis of eigenstates of the oscillator is 

\0„) = I dT{k)0{k) y «(fc>+)_"% -|a(fc,+)|V2 «(fc^"- e-|a(fe,-)|V2|fc)|^^)|^_) (99) 

= m / dT{k)Oa{k,n+,n-)\k,n+,n-) (100) 

where 

Oc«(fc,n+,n_) = ^ n(fc)rv(fc +)"+r.(fe -)"~^~S±l«(fe'±)lV2 (101) 
^Jn+\n-\ 

The average of the 1-oscillator field operator evaluated in such a coherent state is 

{Oa\-Fab{x)\0^)= ydr(fe)e„fa(fe)Z(fe)(a(fc,-)e-^'=-"+^(fe~f)e^'=-^). (102) 

The Poincare transformation of the state implies 

{Oa\Vly-F^b{x)VK,y\Oo:)=K^kb'^{Oo,\-F,d{K-\x-y))\Oo:) (103) 
The coherent-state wave function transforms by 

VK,vOa{k, n+, n_) = e*("++"-)'=-!'e2*("+-"-)®(^''^)0„(A-ifc, n+, n_) (104) 
= O(A-ifc) n e"='=-^e2*'*"=*3(^''=)-^a(A-iA;, s)"=e-l"(^"''^'^)l'/2 (105) 

= O(A-ife) TT -LfTA.a(fc,s)r^e-l^'^.''«('='«)l'/2 (106) 



where 

a(fe, s) ^ TA,ya{k, s) (107) 



is the spin-1 massless unitary representation (|4^). Using this result we get again 

{Oc.\Vly-Fab{x)VK,y\0^) = y"dr(fc)e,b(fc)^(A-ifc)(rA,,«(fc,-)e-*^+7XW^e*^) (108) 

= Ka'Ki\Oo.\-hd{^-\x ~ y))\Oo.) (109) 

showing that those somewhat counter-intuitive forms of ?7a,j, and Va,j^ are consistent with passive Ta^j, transformations 
of classical wave functions. 

With a(/3) and a(/?)^ given by (^)-(|50|) we define the displacement operator 

V{(3) ^ e°('^)''"~°('^) (110) 
= exp(^^ / dr{k){(3{k, s)a{k, s)^ - P{k, s)a{k, s))) (111) 

s •' 

= /rfr(fc)|fc,s)(fc,s|0eS=(^('"''^"'-^^"=) (112) 



which performs a shift of the classical wave function 

Vip)\0^) = |O„+0) (113) 

and commutes with I^'. 

V{I3)^ IkV{l3) = Ik- (114) 
Vacuum states are also coherent states corresponding to a = 0. 

E. Multi-oscillator coherent states 

Consider a family aN{k, s), N = 1,2, . . . of functions and the state 

oo 

= 0\^|Oa„)®.._^|OaJ (115) 
^ N ' 

where 

|0„„)= J dT{k)Oik)\k)\aN{k,+),aN{k,-)} (116) 
and '^'^^iPn = 1- Taking, for example, a^ik, s) — a{k, s)/\/N we find 

{0„\- labium J = I dT{k)eab{k)Z{k)(^a{k,-)e-'''-- + Mk;+)e'''-^) = {Oc\- Fab{x)\Oa) (117) 

i.e. the same result as in the 1-oscillator case. 
A multi-oscillator displacement operator is 

C30 

2(/3) - V{PN)^. . .^V{f3N) = e^^'^)'-^^^ , (118) 

^=1^ N ' 

/3Ar(fe, s) = (3{k, s)/^/N, implying 

^(/3)|OJ = \0^+p) (119) 
P(/3)ta(p, s)P(/3) = a{p, s) + (3(p, s)Lp (120) 

2(/3)UpH/3)-Ip. (121) 



The fact that aN{k,s) = a{k, s) / \/ N will be shown to be of crucial importance for the question of statistics of 
excitations of multi-oscillator coherent states. Let us note that a similar property of coherent states was found in 0| 
when we employed the definition in terms of eigenstates of annihilation operators. 
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F. Multi-oscillator vacua 



Vacuum consists of states with n± = 0, i.e. with all the oscillators in their ground states. Of particular interest, 
due to its simplicity, is the following vacuum state 



\0) = Vm \p) (^. ■ ■(E>\o\ (122) 

where 



N 



\0) = J dr{k)0{k)\k, 0,0) (123) 

Such a vacuum is simultaneously a particular case of a coherent state with a{k, s) = 0. Coherent states are related 
to the vacuum state via the displacement operator 

oo 

2(a) |0) = V^ \Oa^)^. . .<E>\Oc.^\ (124) 
= \OJ. (125) 



G. Normalized 1-photon states 

Consider the vector 

«(/)^IQ). (126) 

Choosing the particular form ( |122| ) we find 

{0\a{f)a{g)^ \0) =Y. I dT{k)Z{k)J{k;s)g{k, s) (127) 
= {fO\gO) {}\g)z. (128) 

fO denotes the pointlike product fO{k,s) = 0{k)f{k,s). Since anyway only the modulus |0(A;)| — Z{k)i occurs in 
the above scalar products one can also work with /s(fe, s) = Z{k)2 f(^k, s). The relation between fs and / resembles 
the one between the bare and renormalized fields jl^. We believe this is more than just an analogy. 
Thinking of bases in the Hilbert space one can take functions fi satisfying 

{Mfj)z = 6,, = {fBAfB,). (129) 



H. Normalization of multi-photon states 



Normalization of multi-photon states is more complicated. In this section we will discuss this point in detail since 
the argument we give is very characteristic for the non-canonical framework. It will be used in Sec. X to show that 
in the thermodynamic limit of a large number of oscillators the non-CCR perturbation theory tends to the CCR one 
but in a version which is automatically regularized. We will also use a similar trick to show that the multi-oscillator 
coherent states have, again in the thermodynamic limit, Poissonian statistics of excitations. 

Denote by the sum over all the permutations of the set {1, . . . , m}. 



Theorem 1. Consider the vacuum state ( |122D with = I for some N. Then 

Jiij^(0|a(/i) . • ■ a(/m)a(5i)^ . . .a(g„)^|0) = ^(/i|g^(i))z ■ . . {frn\ga{ni)) z 



J2 J2 j dT{ki)Z{ki)...dT{k^)Z{k,n)h{ki.si)...f,n{k 
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Proof: The scalar product of two general unnormalized multi-photon states is 
(Q|a(/i)...a(/„)a(5i)^..a(g™)^|0) 

= X! X! / ^r(fci)...dr(fc„)/i(fci,si).../„(fe„,s„)g^(i)(fci,si)...g^(„i)(fcm,sm)(0|/fc^ .../fc^|0) 

cr Sl...Sm 

= ^ ^ / dr(fci) . . .dr(fc„i)/i(A;i, Si) . . . fm{k„i,Sm)ga(l)ikl,Si) . . ■ga{m){km, Sm) 

a S 1 . . . S n^ 

^,{0\...{0\(lk,® + + ...®Ik?j ■■■ [h^ (g>...(g>I+... + I(g>...(g) /fc,„) |0) . . . |0) (130) 



N 



Further analysis of (130) can be simplified by the following notation: 



Nk,=I®...® Ik, 



with j = 1, . . . , m; the sums-integrals 'Y^g. J dr{kj) are denoted by J2k- ■ Then (130) can be written as 



N 



E E /i(fci)---/™(fc"05a(i)(fci)...5.M(fcm)^ E {0\...(0\A,,...Zk^\0)...\0) (131) 
Since m is fixed and we are interested in the limit N ^ oo we can assume that N > m. Each element of the sum 



over ■ ■ ■ in (131) can be associated with a unique point [A^ . . . , Z) in an m-dimensional lattice embedded in 
a cube with edges of length N. 

Of particular interest are those points of the cube, the coordinates of which are all different. Let us denote the 
subset of such points by Cq. For {A, . . . , Z) e Co 

(0| . . . (0| Ak,... Zfc„ |0) . . . \0) = Z{ki) . . . Zikm) (132) 

^ V. ' V ' 

N N 

no matter what N one considers and what are the numerical components in {A, . . . , Z). (This makes sense only for 



N > m; otherwise Cq would be empty). Therefore each element of Cq produces an identical contribution (132) to 
(131). Let us denote the number of points in Cq by A^o- 



The sum can be now written as 



E E f^{^^)---fni{km)9a(l){kl)---ga(m){km)'PoZ{ki)...Z{k,n) 

+ E E /^•••A4M5.(i)(fci)---5.(™)(fc™)^ E {0\...{0\Ak,...Zk^\0)...\0). (133) 

The coefficient T'o = represents a probability of Co in the cube. The elements of the remaining sum over 
(A ... Z) ^ Co can be also grouped into classes according to the values of (0| . . . {0\A}^^ . . . Zk^\0) . . .\0). There are 
TO — 1 such different classes, each class has its associated probability Vj, < j < m— 1, which will appear in the sum 
in an analogous role as Vq. 

The proof is completed by the observation that 



N 



lim -Po = 1, (134) 

I — »-oo 

lim Vj =0, < j. (135) 



Indeed, the probabilities are unchanged if one rescales the cube to [0, 1]"*. The probabilities are computed by means 
of an m-dimensional uniformly distributed measure. N oo corresponds to the continuum limit, and in this limit 
the sets of points of which at least two coordinates are equal are of TO-dimensional measure zero. ■ 
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Comments: (a) T he thermodynamic hmit is naturally equipped with the scalar product yielding orthogonality relation 



of the form (129). However, for small N there will be differences if m is large. On the other hand if N is sufficiently 
large then the values of m for which the corrections are non-negligible must be also large. But then a classical limit 
will be justified and the use of non-canonical coherent states should again give the correct description, (b) Concrete 
values oiV-j for some small m were given in |§. For m = 2: Vq ^ l-l/N ,'Pi ^ l/N] for m = 3: T^o = + 
Vi = 3/N — 3/A^^, 7^2 — l/N'^- In general we do not have to assume that pat = 1. U pn are general probabilities 
then the coefficients involve averages. For m = 2: "Po = 1 - (1/^), Vi = (l/N); for to = 3: Vo^l- (3/7V) + {2/N^), 
Vi = (3/7V) - (3/A^2), V2 = (1/A^^), where (1/A^> = J2nPn/N etc. The normalization in terms of {■\-)z is then 
obtained under the assumption that all those averages vanish, which can hold only approximately, meaning that the 
probability is peaked in a region of large iVs. 

I. States generated by field operators and a first hint indicating absence of ultraviolet divergences 

Consider the single-oscillator vector potential operator Aa{x) which is related to ( |68|) by A^{x) — ®^Aa{x) and 
acting with Aa{x) on a single-oscillator vacuum \0) define the vector 

\Aaix)) ^ Aa{x)\0) (136) 
dr(fe)e*^ (ma{k)a{k, -)^ + ma{k)a{k, +)^) \0) (137) 

= -i j dT{k)e'''-''OM{k)\k)(n2aik)al\0,0) +ma{k)al\0,0)^ (138) 

Its multi-particle analogue is 

\Aai^)) ^ Aai^m) (139) 

= (B]^^,^{\Aa{x))\0)_^+. . .+ \0)_^\Aa{x))). (140) 

N-1 N-1 

The positive definite scalar product 

{Aa{y)\i~9'"'Mt{^)) = {A^m-g'^'M^ix)) = 2 / dr(fc)e*(--^)z(fc) (wi) 



shows that there is no ultraviolet divergence at a; = y since / dr{k)Z{k) = 1. 

It is easy to understand that the same property will hold also for general states. To see this let us write the 
single-oscillator field operator as a function of the operator ka = J dT{k)ka\k){k\, i.e. 

Aaix) = iTOQ(fc)(e~''''^ (g) a+ - e''''^ «) al) + ima{k){e^ ''''''' ® - e''''^ ® a\). (142) 

The operators ma{k) and ma{k) are functions of the operator k and are defined in the standard way via the spectral 
theorem. Moreover, they are complex combinations of two bounded operators representing directions of transverse 
polarizations. The remaining operators {aj, a|, and e*'*^'^) are also well behaved. Particularly striking is the fact 
that the distribution J dr(fe)e*'^'^ is replaced by the unitary operator 



dr(fc)e*^|fc)(fc|. (143) 



The latter property is at the very heart of the regularities encountered in the non-canonical formalism. 

Now, it is widely known that configuration-space renormalization of ultraviolet divergences can be re duced to an 



appropriate treatment of products of field operators on the diagonals x = y |16yi71| . The formula (141) is a strong 
indication that such divergences may be absent in the non-canonical framework. Actually, the analysis of perturbation 
theory in nonrelativistic quantum optics given in [|l| and further elaborated in Sec. X of the present paper, shows 
that Z{k) occurs in exactly those places where ultraviolet form-factors are expected to appear. The same property 
will hold for non-canonically quantized fermionic fields. 
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VII. STATISTICS OF EXCITATIONS 



It is an experimental fact that laser beams produce Poissonian statistics of photocounts. At the theoretical level of 
canonical quantum optics the Poisson distribution follows trivially from the form of canonical coherent states. In the 
non-canonical case the exact Poisson statistics is characteristic of the single-oscillator {N = 1) sector. For 1 < iV < oo 
the statistics of excitations is non-Poissonian. At the other extreme is the thermodynamic limit for multi-oscillator 
states. In what follows we will show that in the limit N oo one recovers the same Poisson distribution as for A'^ = 1. 
This, at a first glance unexpected, result justifying our definitions in terms of displacement operators is a consequence 
of certain classical universality properties of the Poisson distribution. 

We will also return to the question of thermal states and the Planck formula. In ^ it was argued that non-CCR 
quantization implies deviations from the black-body law. However, a consistent interepretation of the field in terms 
of the thermodynamic limit shows that no deviations should be expected. 



A. Multi-oscillator coherent states 



To study the thermodynamic limit of multi-oscillator coherent states we simplify the discussion by taking an exactly 



A^-oscillator coherent state (115) (A^ ^ 1 is fixed and = 1), i.e. 

|OJ = |Oa„)0...0|Oo„) (144) 



JV 



where 



|0„„)= J dT{k)0{k)\k,s)\a{k,+)/VN)\a{k,-)/VN). (145) 
The average number of excitations in this state is 

{n)=Y, j dT{k)Z{k)\a{k,s)\\ (146) 

The simplest case is the one where a{k,s) ^ a = const. Then (n) — |ap and the statistics of excitations of 
single-oscillator coherent states |Oq„) is Poissonian with the distribution p„ = g"'""' \aM\^"' /n\ 

m excitations distributed in the ensemble of A^ oscillators can be represented by the ordered m-tuple (ji, . . . ,jm), 
I < ji < ■ ■ ■ < jm < N. For example, for m = 10, A^ = 12, the point (2,2,2,5,5,7,7,7,11,11) represents 10 
excitations distributed in the ensemble of 12 oscillators as follows: 3 excitations in 2nd oscillator, 2 in the 5th one, 
3 in the 7th, and 2 in the 11th. Such points form a subset of the cube [0, A^]™, the interior of the set corresponding 
to points whose all the indices are different. The latter means that the interior represents situations where there 
are m oscillators excited, and each of them is in the first excited state. The boundary of this set consists of points 
representing at least one oscillator in a higher excited state. Probabilities of events represented by points with the 
same numbers of repeated indices must be identical due to symmetries. Intuitively, the Poissonian statistics of the 
thermodynamic limit follows from the fact that the probability of finding a point belonging to the boundary tends 
to zero as A^ increases. The statistics is dominated by Bernoulli-type processes with probabilities related to the two 
lowest energy levels of a single oscillator in a coherent state. 

To make the argument more formal we introduce the following notation: 

aW = {x e iV";TO > l,x^ (ji, . . . , j™), 1 < < . . . < < TV} 
^iV-n^ = {a; e A,(f);a; = (n, . . . , ii, . . . i^, . . . ,ife),n < . . . < ik} 

U (147) 

(rii...nfc)5^(l---l) 

If we add a single-element set Ag^-* containing the event representing A^ oscillators in their ground states we can 
represent the set of all the events by the disjoint sum 

CO 

A(^) = U Xi^) (148) 

m=0 
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The probability of finding the partition m = ni + . . . + rife is 



ni...nk I I 

ni!...nfe! 



where Nni...nk is the number of elements of X^^.n^ C The probability that m excitations are found is 



(149) 



(150) 



ni...nk 



the sum being over all the partitions of m. Denote by P{Yji^^\Xm^) the conditional probability of finding at least 
one oscillator in the 2nd or higher excited state under the condition that the sum of excitations is m > 1. We first 
prove the following 



Lemma 1. 



hm p(yW|xW)=o. 

JV— >CX) 



Proof: Since vi^^ n X^^ = vi^^ one finds 



) P{Xnl}.n^-) 



(Ttl...7tfc)7^(l...l 



< 



(ni...nfc)yi(l...l) 

l + iVl...l( E 

(ni...nfc)^(l...l) 



However, 



I](ni...Kfc)#(l...l) ^ni...nk ^ ^ 

on the basis of the geometric argument we gave in the previous section. This ends the proof. 
The main result of this section is the following version of the well known Poisson theorem: 

Theorem 2. Assume that a{k, s) = a = const. Then 



(151) 



(152) 



(153) 



lim P(X(f)) 

AT— >cx) 



e-l"l lap™ 



Proof: As an immediate consequence of the lemma we find 



lim P(X(f )) = hm P(X(f ) - r^)) 

N^oo N^oo 



(154) 



(155) 



which means that in the thermodynamic limit we can treat excitations of the oscillators to 2nd and higher excited levels 
as events whose probability is zero. The probabilities of ground and first excited states follow from the single-oscillator 
Poisson distributions but conditioned by the fact that only the lowest two levels are taken into account. 
We thus arrive at the standard Poisson process with 



Pn = 



Pi 



\aN\ 



H'/N 

PQ+Pi l+|ajv|2 l + \a\yN 



(156) 



and limjv^cx) ^Pn = laP. 
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B. Thermal states 



A single-oscillator free-field Hamiltonian H has the usual eigenvalues 

E{LU,n) ^hu;(n + ^y (157) 

The eigenvalues of the free-field Hamiltonian at the multi-oscillator level are sums of the single-oscillator ones. In 
101 it was assumed that the Boltzmann-Gibbs distribution of thermal radiation should be constructed in terms of I£. 
Let us note, however, that such a construction makes use of as if it was a Hamiltonian of a single element of a 
statistical ensemble. The discussion of the thermodynamic limit we have given above, as well as the results of 
suggest that H_ is the Hamiltonian of the entire ensemble of systems described by H, and it is H and not which 
should be used in the Boltzmann-Gibbs distribution. Then, of course, the result will be the standard one and no 
deviations from the Planck formula will occur. 

VIII. COMMUTATORS OF 4-POTENTIALS AND LOCALITY 

A straightforward calculation shows that the multi-oscillator vector potential operator satisfies, for any space-time 
points X, y, z, the commutators 



+ y"(ir(fc)/fe(m„(fc)77Zfc(fc)e*'=-(^-^) -ma(fe)mb(fc)e*'=-(^"'')) (158) 



To obtain more in sigh t as to the meaning of the commutator (15S) consider its coherent-state average evaluated in a 
state of the form ( |115| ): 

{OJ[Aaix),A,{y)]\0^) = I dT{k)Z{k)(ma{k)fh,{k)e'^<y--^ - m,(fc)™f,(fe)e*(^-^)) (160) 

dT{k)Z{k) (m„(fc)mfc(fc)e*'=-(^-^) - ma(A;)77i6(fc)e*'=-(^~y)) . (161) 

The Minkowski-space metric tensor can be decomposed Q in terms of null tetrads as follows 

gab ^ K(^b + ^^ah - rnafhb ~ rrihffia. (162) 
With the help of this identity we can write 

(OJ[AJa;),^(y)]|OJ = J dT{k)Zik){kaMk) + kbCOa{k)){e^'<y-^^-^'<^-y^) 

+gab f dr(fc)Z(fc)(e'^-("-^) - e''^-^^-")) (163) 



It is known that terms such as the first integral vanish if the potential couples to a conserved current. The same 
property guarantees gauge independence of the formalism. Therefore we can concentrate only on the explicitly gauge 
independent term proportional to gab- Denote 

Dz{x) = I J fir(fc)Z(fc)(e-*^ - e*'=-^) (164) 

For Z{k) = const = Z we get Dz{x) proportional to the Jordan-Pauli hmction, 

Dz{x)^ZD{x) (165) 

which vanishes for spacelike x. However, the choice of constant Zik) is excluded by the requirement of square- 
integrability of O. Therefore the requirement that the vacuum state be square-integrable seems to introduce some 
kind of nonlocality into the formalism. 
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There are two possibilities one can contemplate. First of all, one can perform the calculations with arbitrary O 
and then perform a renormalization (we have seen that such a step is necessary even in the nonrelativistic case). 
After the renormalization we can go to the "flat" pointwise limit Z{k) — ^ 0, || O ||= 1, corresponding to the uniform 
distribution of all the frequencies. Second, performing the calculations in a preferred frame we can consider equal-time 
commutation relations 

(OJ[AJt,a:),^(i,y)]|OJ = j dr(fc)Z(fe)(fc„..,(fc) + fc,c.,(fc))(e*(^^--T-e*(-~-^~)) 

+9ab j dr(fc)Z(fc)(e*(^-^^) - e'^-(«^-^)) (166) 

The last term will vanish if 

Z{k) = Z{~k) (167) 

i.e. if the vacuum is 3-inversion invariant. This can hold, however only in one reference frame unless O is constant, 
which we exclude. 

One can conclude that non-canonically quantized electrodynamics is not a local quantum field theory, at least in 
the strict standard sense. This is not very surprizing since the presence of Z{k) in the integrals introduces some kinds 
of effective extended structures, a consequence of nontrivial structures of non-canonical vacua. The issue requires 
further studies. In particular, it is important to understand an influence of the thermodynamic limit N oo on 
locality problems in the context of relativistic perturbation theory. 

There is an intriguing analogy between the kind of non-locality we have obtained and the one encountered in field 
theory in non-commutative space-time [0-01. 



IX. RADIATION FIELDS ASSOCIATED WITH CLASSICAL CURRENTS 



The problem of radiation fields is interesting for several reasons. First of all, the radiation fields satisfy homogeneous 
Maxwell equations so that the theory we have developed can be directly applied. Second, this is one of the simplest 
ways of addressing the question of infrared divergences within the non-canonical framework. 

It is widely known [|2 19 1 that in the canonical theory the scattering matrix corresponding to radiation fields 
produced by a classical transversal current is given, up to a phase, by a coherent-state displacement operator 
yJ{v)-A,n{y) ^ Qj^g q£ consequences of such an approach is the Poissonian statistics of photons emitted 
by classical currents. An unwanted by-product of the construction is the infrared catastrophe. 

Starting with the same S matrix but expressed in terms of non-canonical "in" fields we obtain the non-canonical 
displacement operator. Photon statistics is Poissonian in the thermodynamic limit but the infrared catastrophe is 
automatically eliminated. 



A. Classical radiation field 

Let us assume that we deal with a classical transversal current Ja(x) whose Fourier transform is Ja{k) — 
J d'^xe^^'^ Ja{x). Transversality means here that 

J,(|fc|,fc) =m„(fc)Jw(|fc|,fc) + m,(fc)Joi,(|fc|,fe). (168) 

Formally, a solution of Maxwell equations 

UA^{,x) = Ja{x)L (169) 

can be written as 

Mx) = A^,^{x) + j d^yD,,tix - y)Ja{y)L (170) 

= Aaontix) + / d^VD.A.{x - y)Ja{y)L (171) 



Here Aain ^'^'^ Aaout s-^^ solutions of homogeneous equations. D^ct and I?adv are the retarded and advanced Green 
functions whose difference is the Jordan-Pauli function 
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D{x) = i j dr(fc)(e-''^"^ - e*^). (172) 
The 4-potential of the radiation field is 

Aar.di^) = Aaouti^) " Aaini^) = J d^Dix - y)Ja{y)L (173) 

and leads to the field spinors 

£x'y'rad(^) = / rfr(fc)(^(;f,|^^Jy|y,)(fc)e*- +^(X'|^^Jy|yo(fc)e^*')l- (175) 
Comparing these formulas with expressions (^ and (|^) valid for all solutions of free Maxwell equations one finds 

f{k,+) = ~m'^{k)Ja{\k\,k) = Jov{\k\,k) = j{k,+) (176) 

/(fe, -) = -m%k)Ja{\klk) - Jio'(|fe|, fe) = j(fc, -) (177) 

a(A;, s)out = a{k, s)in + j(fc, s)L (178) 



B. Non-canonical radiation field 



Formula (|l78|) is analogous to the one from the canonical theory. It is clear that although a(fe, s)in and a(fc, s) 



out 



cannot be sim ultan eously of the form given by (53) and (|58|), they do satisfy the non-CCR algebra (|65|). In spite of 
this the result ( 178 ) is not very satisfactory. Indeed, one expects that the scattering matrix describing a non-canonical 
quantum field interacting with a classical current is 

5 = e^*e"*/''''''^(^^-^'"^^) (179) 

with some phase </>. Then 

AaoA^) - S^Aain{^)S (180) 

= Aainix) - * / d^Vj\y)[Aain{x),AbM- (181) 



Employing (|15§), ( |176D , ( p7^ ) one can write 
where ma = ma{k), and 

a(fc, s)out = a{k, s)in + j(A;, s)Lk (184) 
= P(jf a(fc,s)i„2(j). (185) 

Consequently, the S matrix is in the non-canonical theory proportional to the non-canonical displacement operator 

S^e"l"D{]). (186) 

This fact will be shown to eliminate the infrared catastrophe. 
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C. Propagators 



Evaluating the average of (181) in a coherent state \0_^) one finds 



(OaUarad(a^)IO J = J d%Dz{x - y)Ja{v) + gauge term. (187) 



The irrelevant gauge term is a remainder of the first part of (163). As expected the radiation field does not depend 
on what a one takes in the coherent state, but does depend on the vacuum structure. The presence of the regularized 
function Dz{x — y) instead of D(x — y) implies that the radiation signal propagates in a neighborhood of the light cone. 
Any deviation from c in velocity of signal propagation can be regarded as an indication of a non-constant vacuum 
wave function 0{k). To have a feel of the scale of the nonlocality assume that 0{k) is constant up to, roughly, the 
Planck scale. Corrections of the order of the classical electron radius can be seen in the distance travelled by light 
if the photon travels 100 light years. Moreover, even in the orthodox canonical quantum electrodynamics a detailed 
analysis of signal propagation leads to small deviations from velocity of light, especially at small distances It 
may be difficult to experimentally distinguish between the two effects. A similar effect was predicted for Maxwell 
fields in non-c ommutative spac e-tiin e [^,0]. 
Using ( |162| ) one can rewrite (|l58| ) as 



[Aai^),Abiy)] - 9ab J dnk)I_^(e^^<--y) - e^^<y-^)^ +... (188) 

where the dots stand for all the terms which are gauge dependent and do not contribute to physically meaningful 
quantities. We can therefore identify 

D{x) = i J dT{k)lf,(e-'^--' - e'^-^^ (189) 
as the operator responsible for the appearence of the smeared out Jordan-Pauli function Dz in the coherent-state 



average (187). Djx) is a translation-invariant scalar-field operator solution of the d'Alembert equation, i.e. 



nD{x) = 0, (190) 
V\^yD{x)V^^y^D{K-'x). (191) 



The operator analogues of retarded and advanced Green functions are 



D,,,{x) = Q{xo)D{x), (192) 
^advC^) = -Qi-^oM^), (193) 
D{^) = D,.,ti^) - D^^^ix). (194) 



Eq. (19C) implies that the operators 



I d^yD^^,{x - y)Ja{y) (195) 

and 

J d^yD^^^ix - y)Jaiy) (196) 

differ by at most a solution of the homogeneous equation 

UA^ix) = 0. (197) 



( 190 ) implies also that one can define 

m ~ DD,^^{x) = aD,,,{x). (198) 



It follows that having a solution ^^[^(a;) of (197) one can define another solution A^^^^^{x) of (197) by means of 
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Aa[x) = A^,^{x) + j d^yD^^,[x ~ y)Ja{y) (199) 
= A,„,t(x) + / d^yD,^,{x - y)Jaiy), (200) 



simulatneously guaranteeing that the correct 5- matrix conditions (179), ( |l80| ) are fulfilled up to, perhaps, a gauge 
transformation. A^{x) is a solution of 



where 



aA,{x) = J,{x) (201) 



JJx) = / d''ySix-y)Jaiy). (202) 



D. The problem of infrared catastrophe 

To close this part of the discussion let us consider the issue of infrared catastrophe. We have to compute an avarage 
number of photons in the state D{j)\Q}- The number-of-photons operator is 

n = ®i(l®^ajas). (203) 

S 

The 1 in the above formula is the identity in the k space and as satisfy CCR. The average reads 

(n) = (0,N0,) =E / dmZ{k)\jik,sr. (204) 
The four-momentum of the radiation field is 

(Pa) - (0,|PJ0,> - ^ / dr(fc)fc„Z(fc)|j(fe,s)|2. (205) 

0{k) belongs to a carrier space of an appropriate unitary representation of the Poincare group. As such this is a 
differentiable function vanishing at the origin fc = of the light cone. This is a consequence of the fact that the 
cases fc = and A: 7^ 0, A:^ = 0, correspond to representations of the Poincare group induced from SL{2, C) and -E(2), 
respectively (for another justification see |Tl|). 

Hence, the regularization of the infrared divergence is implied by relativistic properties of the field. It is quite 
remarkable that all the divergences are regularized automatically by the same property of the formalism: The nontrivial 
structure of the vacuum state. In the case of ultraviolet and vacuum divergences the regularization is a consequence 
of square integrability of O. 

X. SELECTED QUESTIONS OF PHYSICAL INTERPRETATION 

In this section we want to address three problems which allow to better understand on physical grounds the formal 
structures characteristic of our choice of field operators. First we will discuss in a simplified model (scalar field, 
discrete spectrum of frequencies) the relationship between the non-CCR approach and the canonical one. We will 
perform this in two versions: One using an exact result, and the second one based on perturbation theory. Making 
use of the same model we will interpret the factor I/^/N occuring in the definition of multi-oscillator field operators, 
and finally we will look more closely at the structure of the state space. 
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A. Canonical limit and the second hint indicating absence of ultraviolet divergences: Exact example 



We now want to make a link between the non-CCR theory with nonunique vacuum and the one based on CCR and 
the unique vacuum. One may regard it as a test of experimental consequences of the non-canonical approach. The 
analysis of Poisson statistics of coherent states indicates that the link may be provided by the thermodynamic limit 

— !■ cx). A similar conclusion follows from the analysis of orthogonality properties of multi-photon states. 

In what follows we shall discuss this in more detail on a simplified model of a non-CCR scalar field in dipole 

approximation interacting with a two-level atom located at a; = 0. It is known that without a regularization of the 

1 1 

coupling parameters by gk 9kZ^, with for |fe| oo, the model is highly ultraviolet divergent (even more 

than its relativistic counterpart). On the other hand the vast majority of experimental tests of quantum radiation 
fields is based on this type of calculation. 

We will compare predictions of two models: The CCR one with regularized coupling constants, and the non-CCR 
quantized field without regularization of the Hamiltonian. 

Consider the interaction-picture Hamiltonian 



H{t) - E V 2^ K"'^^"'^'* - R-a^^e^^"') . (206) 



. ha I h 

The coupling constant a depends linearly on the ratio eo/^/V (charge eq and quantization volume V). We assume 
the CCR algebra [afc,aj[,] = 5kk' with the unique vacuum |Occr), and R± = |±)(=f|. 
Of particular interest is the vacuum-to- vacuum amplitude 

Fit) ^ {0\Uit)\0) (207) 

which can be exactly computed in both models. In the CCR case, taking |0) = |Occr)|+), in order to compute the 
amplitude the Hamiltonian H{t) has to be replaced by its regularized version 



Evaluating the vacuum average of the Dyson expansion one finds the Volterra-type integral equation 

Fit) = 1-C f dt^ r dt^fzih - t2)F{t2) (209) 



where C = a'^h/8 and 



fzir) = y2Zk . (210) 



The amplitude expressed in terms of the Laplace transform, with integration along a contour located to the right of 
all the poles of the integrand, is 

Fit) = -L /■ dz , —. (211) 

Let us now take the interaction-picture Hamiltonian 



k 

,t 1 _ 



mt) = E J^i^+^ke-'^"' - R-ak^e^^"') , (212) 



2uik 

with the non-CCR operators [a/,, a^,] = Skk'Lky "^ki-k = L, and the non-canonical vacuum 



\0)^^^/m E Ofei---Ofe.Jfci...fejv) (213) 

N=l fci...fcjv 



which is of the form (122). Let us note that now we do not introduce any regularization of the coupling parameters. 
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Setting Z-k = |Ofcp and denoting the non-canonical amplitude by F' one can write 



F'it)=J2PN Zk,...Zk^{ki...kM\U{t)\ki...kN) 

N^l ki...kN 

oo 

= 'Ypn ^ Zfe^ . . . Zfe„F(i)fc^...fc„. 

N=l fei...fcjv 

The Dyson expansion leads to a similar Volterra-type equation as before 



N In 



dt2f{tl-t2)k,...k«F{t2) fei ...fcjv 



with 



f{T)ki...kN - X! 

LOkj — \kj\, and the same constant C. The solution is 



e 3 



(214) 
(215) 

(216) 
(217) 

(218) 



where F is any c onto ur parallel to the imaginary axis and to the right of all the poles of the integrand. 

The solution (218) looks very sim ilar to ( ^ll[) but there are also evident differences. For example the sum in the 
denominator of the inegrand in (21S) is finite, involves no regularizations, and C is replaced by C/N . 

To see the link between F{t) and F'{t) let us assume that there are exactly N oscillators (i.e. pn = 1 for some 
N and zero otherwise) and consider the thermodynamic limit N oo. This is the same type of reasoning we have 
employed in discussion of the Poisson statistics of coherent states and normalization of multi-photon states. 



Theorem 3. Under the above assumptions 



lim F'{t) = F{t). 



(219) 



Proof: The poles of the integrand in (218) are equal to the eigenvalues of 

... 



c 



C 
ojjvJV 



iA,. 











(220) 



and, hence, are purely imaginary. The parameters C, TV and toi, . . .lon in (218) are fixed, integration is over any 
contour localized to the right of all the poles, and the poles are imaginary. We assume that the spectrum of ojs 
contains a minimal u; > (as in cavity) . It follows that the contour can be shifted sufficiently far to the right so that 
the inequality 



1 



N 

\N ^ ujj z i A^^- + z 



is satisfied for any choice of , . . . , w^v , and 



< 1 



N 



(221) 



N l-'j = \ ujj z iA.^.+z 71=0 j = l J 



(222) 



The amplitude of interest can be thus written as 
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= ^ E ^'^-•■^fc"E(-^)" / ^^4TT^(E--7r^) (223) 

ki...kN n=o j = i ^ 

1 ^ /■ ^* 1 ^ 1 1 



+ 2^ E ^fei---^fc«^,;>ivWfei...fc«- (224) 

The convergence of the geometric series guarantees that for any e there exists such that 

{t)k,...k.\ < e. (225) 

fcl...fciv. 



Therefore 

Urn F'(i) = — hm y(-C)" / V Z^, ...Zfe„ — fV 

ri=0 ^ ki...kN j = l ^ 



(226) 



f%Vo + ■■■ + ffiVn] , (228) 



hm / dz — 



27ri N^oo ./r 



where Vq . . .Vn are the probabihties employed in the proof of Theorem 1 and Hmjv^oo < oo for aU j. Using 
hnijv^oo T'o = 1 and hm^r^oo Vj = for j > we find 

lim F'U) = — f dz — ^ — 1 . (229) 



It is important to note that F{t) obtained as the Umiting value of F'{t) contains Zk = |OfcP satisfying the normal- 
ization condition J2k = 1- Therefore although the thermodynamic limit regularizes the non-canonical amplitude 
as if it originated from the regularized canonical Hamiltonian i?reg) the regularization involves functions Z^ which are 
somewhat unusual from the viewpoint of standard quantum optics. 

Indeed, the typical choice one finds in the literature is of the form 

^ ' + A2 \ for |fe| > A ^ ^ 

and it is essential for the standard interpretation that Zk ^ 1 with A oo. The regularization we arrive at must be 
square summable (or square integrable) and normalized. So assume 

Z,(Ai, A2) = ( f = 7"? \t\ I Sa^'M (231) 
^' \ tendstoO for |fc| ^ [Ai,A2] ^ ' 

But since Zk — 1, the limit A2 — > 00 is accompanied by Z ^ 0. An agreement with the standard formalism is then 
obtained if one performs the latter limit under the constraint CZ = Ccxp, i.e. renormalizes the coupling parameter a 
by ^exp — cxZ 2 ^ which is essentially the charge renormalization Ccxp — ^qZ'^ . 
All the examples studied so far suggest the following 

Conjecture: The thermodynamic limit N —i- 00 maps predictions of the non-canonical theory into those of the 
canonical theory regularized by square-integrable formfactors. Elimination of regularization must be accompanied by 
charge renormalization. 



23 



B. Canonical limit and the third hint suggesting absence of ultraviolet divergences: Perturbation theory 

Let us now discuss the thermodynamic Umit of a perturbative expansion to order n of an arbitrary ampUtude 

FMt) - (*/|C/(t)|*,). (232) 

Examples of low order perturbation theory in concrete examples were explicitly trea ted in . Here we want to discuss 
an arbitrary order of a general amplitude corresponding to the Hamiltonian ( ^12[ ) which has led to the Conjecture 
from the previous section. It is sufhcient to discuss the amplitude for the general multiphoton states discussed in 
Theorem 1. 

The interaction Hamiltonian ( ^12[ ) is applicable to fields with discrete frequency (i.e. cavity) spectrum. This 
is technically a simple case since momentum eigenvectors are normalized in volume V by the Kroenecker deltas. 
Performing the expansion to order n of the amplitude we arrive, for each Feynman diagram in momentum space, at 
the term proportional to 

(Q|a(/i) • • .a(/,„)Mono(a(fci), . . . , a(fcn))a(gi)^ . . . a(gm')^|0) (233) 

where Mono(. . .) is a m onon omial of order n in creation and annihilation operators. Due to the Kroenecker-delta 
normalization the term ( ^33| ) vanishes if and only if an analogous term does in the canonical theory (cf. the discusion 
in 0). For this reason the non-CCR amplitude produces the same Feynman diagrams, but the Feynman rules are 
modified by the non-CCR algebra. 

The latter modificatios are easy to identify. Simply, each time a commutator of creation and annihilation operators 
produces a delta in the standard theory, say Skk', here it is replaced by Skk'Lk- The Wick-theorem- type manipulations 
one has to perform under the vacuum average are of the same type as those we have met in the proof of Theorem 1 
while normalizing multi-photon states. 

The essential and final step follows from the main property of the thermodynamic limit: limjv^ooT'o = 1 E^nd 
limjv— ►oo "Pj = for j > 0. Therefore only those terms count in the limit which have each vacuum average of the 
product . . . /j,^ replaced by Zk^ . . . Zk„ ■ Choosing the initial and final wavepackets fj , gj normalized with respect 
to the scalar product {■\-)z from Theorem 1 we find that the modificaltion of the amplitude is of precisely the type we 

have formulated in the Conjecture: Coupling constants in momentum space are regularized according to gk — > gkZ^ 
and the ultraviolet infinities are regularized by square summability of Zk- 



C. The physical meaning of the parameters occuring in the multi-oscillator non-CCR representation 



The Hamiltonian ( pi 2]) can be also written as 

H{t) =^^y] \f^(R+ a^e-'"^-' - i?_ ®^ a^^e^^'^*) . (234) 

h 

The choice oi aN = was motivated by purely formal reasons. With this choice: (1) the RHS of field commutators 
satisfied the resolution of unity, (2) the statistics of coherent-state excitations was Poissonian in the thermodynamic 
limit, (3) the formfactors Z{k) were regularizing the interaction term in the thermodynamic limit, (4) non-canonical 
coherent-state averages at the multi-oscillator and single-oscillator levels were identical, (5) normalization of single- 
photon and multi-photon states could be given in terms of the same scalar product {•\-)z in the thermodynamic 
limit. 

The restriction H_{t)N of H_{t) to an A^-oscillator subspace reads 



N 



j = l ^ k 

(i) 

where — I(^...ak---'^I with the single-oscillator annihilation operator at the jth place. 

It is clear that H_{t) jv represents the interaction Hamiltonian of a single 2-level system interacting with N indefinite- 
frequency harmonic oscillators. The factor means that the strength of the coupling decreases with growing N, 
a property whose physical interpretation is well known and occurs, for example, in the Hepp-Lieb treatment of the 
Dicke model [||j23|. 
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The Dicke model is in a sense dual to the one we consider. It represents a single harmonic oscillator interacting 
with N 2-level systems via the interaction-picture Hamiltonian 



H[t)NBicke - 2_^i-7^^y -^yn+ ane - H_ a^^e 



(236) 

^/ ;v z V zsi V ■ / 

Here r'£' = / . . . i?± . . . /, and [aj^, ajj = /. 

Now, assuming the particular representation of \a n , a|^] = / given hy fl = Ci;fc|fc)(fe|® 1, / = 1® 1, [ofe, aj] = 5kil, 
o-k = \k){k\ ® cik, an = J^k '^k, we can rewrite (236) as 



N 



J— 1 fc 



Although the Hamiltonian (236) is formally a single-mode one, the form ( [237] ) resembles the non-canonical interaction 
term discussed in the previous subsection. What is important, the TV-dependent coupling constant can be derived on 
physical grounds if one assumes a constant density N/V of (bosonic) atoms. To see this, one writes the parameters 
of the Dicke-Hepp-Lieb Hamiltonian as 

a{h/2f/'^ = ujady/2T:hp (238) 

where p — N/V is the density of atoms, d the dipole moment, and luq the atomic frequency. The terniodynamic limit 
of the Hepp-Lieb approach is performed under the assumption a ~ const which is equivalent to the constant density 
P- 

This is precisely what happens in our approach to the thermodynamic limit. It follows that the choice of the 
non-CCR algebra whose RHS satisfies the resolution of unity is physically equivalent to the requirement that the 
electromagnetic field consists of indefinite-frequency oscillators of constant density N/V. The thermodynamic limit 
is then equivalent to the infinite volume limit V ^ oo. 

Finally, let us note that taking the direct sum over N we find 



WDickc 



where ^_|_ = (B^R±, i.e. a structure analogous to the interaction Hamiltonian in non-canonical quantum optics. 



D. Structure of the space of states 



The Conjecture we have formulated above implies that probabilities have to be computed as if the vacuum was 
unique. However, it is essential for the non-canonical construction to have an infinite dimensional space of different 
vacua. Otherwise it would be impossible to associate the spectrum of frequencies with a single harmonic oscillator. 
There is no contradiction if we treat the Fock space genereated from a vacuum state \0) as a fiber over \0) in a vector 
bundle with Fock fibers and the base space consisting of vacua. Vacuum fiuctuations represented by the zero-energy 
part of the Hamiltonian generate a motion (a flow) in the base space and thus play a role of bundle connection. The 
"vacuum picture" corresponds to a choice of connection. 



XI. CONCLUSIONS AND FURTHER PERSPECTIVES 



Fields with the property of having at the RHS of the commutator an element from a nontrivial center of the algebra 
are usually termed the generalized free fields. In the context of our formalism we find this term misleading for two 
main reasons. First of all the generalized free fields are "free" whereas the examples we discussed include scattering of 
radiation, spontaneous emission and interaction picture perturbation theory. All these processes involve interactions 
with charges. Secondly, we arrived at the particular non-CCR representation on the basis of a concrete physical model 
of the indefinite frequency oscillator. The mathematical structure of the model does not follow the typical limitations 
imposed on generalized free fields. For example, the vacuum states are Poincare non-invariant. The same concerns 
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the RHS of the non-CCR algebra which is only Poincare covariant and satisfies the resolution of unity. The issues of 
locality have to be formulated in terms of certain operator generalizations of Green and Jordan-Pauli functions which 
lead to very special locality properties and eliminate problems with multiplication of field operators at the same point 
in configuration space. The notion of the thermodynamic limit is typical of the representation we use and does not 
seem to have been investigated in the context of axiomatic quantum fields. 

The construction we are developing has a slightly different logic than the usual axiomatic quantum field theory. It 
seems that what we are doing is somehow in-between the generalized free fields, nonlocal quantum field theory, and 
noncommutative geometry. 

An issue which has not been addressed so far is how to quantize fermions. Some results on the Dirac equation are 
known already and will be presented in a separate paper. 

Another problem is to embed the concrete non-canonical quantization procedure we have proposed into a more 
abstract scheme of quantizations in a C* -algebraic setting. The fact that the right-hand-side of commutation relations 
is not an identity but rather an operator belonging to the center of the algebra suggests directions for generalizations. 
It seems there is a link with the work of Streater on non-abelian cocycles [|4j . An appropriate version of a coherent- 
state quantization based on the formalism of Naudts and Kuna ||25|l is in preparation. 
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